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Abstract 

Di Paola and Falsonc's formula is widely used in studying stochastic 
dynamics of nonlinear systems under Poisson white noise. In this short 
communication, an alternative expression is presented. Compared to 
Di Paola and Falsonc's original expression, the alternative one is appli- 
cable under more general condition, and shows significantly improved 
performance in numerical implementation. The alternative expression 
turns out to be a special case of the Marcus integrals. 

Keywords: Stochastic dynamics; Poisson white noise; Non-Gaussian 
white noise; Nonlinear systems; random vibration. 

1 Introduction 

In the presence of random effects, dynamics of systems in engineering and 
science is often modeled by stochastic differential equatinos (SDEs) driven by 
Gaussian or non-Gaussian processes. The solutions of SDEs are often denned 
by Ito or Stratonovich integrals. When the excitation is parametric Gaussian 
white noise, Stratonovich SDEs are different from the corresponding Ito 
SDEs: they have an additional term, the so-called WZ correction term [13\ 
[H] . The studies in [7J [15] attribute the extra term to the conversion from 
physical white noise to ideal white noise. For Poisson white noise excitation, 
the correction terms for converting Stratonovich SDEs into Ito SDEs are 
first presented by Di Paola and Falsone [U 0] • Although some disagreement 
has been brought up concerning Di Paola and Falsone's correction term in 
earlier days [U[5], the Stratonovich SDEs proposed in [31 0] have been verified 
extensively and are widely used to study stochastic dynamical systems under 
Poisson white noise excitation [161 [T71 [T2"l [2] . 

The main objective of this short communication is to give an alternative 
expression of Di Paola and Falsone's formula. Compared with the original 
one, the alternative formula is applicable under more general condition and 
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shows better performance in numerical computation. Note that we con- 
sider only the one-dimensional cases. The conclusion can be generalized to 
high-dimensional cases which, however, will not be discussed in this short 
communication. 



2 Review of Di Paola and Falsone's formula 

Consider a dynamical system governed by 

Z(t) = f(Z(t),t)+g(Z(t),t)C(t), (1) 

where Z represents the state variable, / and g are deterministic functions 
of Z and t, C(t) is Poisson white noise expressed as 

N(t) 

C(t) = Y,RkS(t-t k ). (2) 
fc=l 

In Eq. ([2]), N(t) is a Poisson process with intensity parameter A, S(t — tk) is a 
Dirac function at tf., is a random variable representing the fe-th impulse. 
The Poisson white noise C(t) is regarded as the formal derivative of some 
compound Poisson process C(i), 

N(t) 

C(t) = Y,RkU(t-t k ), (3) 

fc=l 

where U is a the step function. 

In Ito sense, ([1]) is often written in the form 

dZ(t) = f(Z(t),t) dt + g(Z(t),t) dC(t), (4) 

Di Paola and Falsone suggest to interpret di]) as [31 S] 
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dZ(t) = f(Z(t-),t) dt + ^ g0)(Z ff — (dC(t)y , (5) 

where 

' gW(x,t)=g(x,t), 

g^(x,t)=g(x,t)^- 1 \x,t), fori > 2 



i=i J 



(6) 
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© is often called the Stratonovich SDE for © [3j H EE HZ]- Note that 
dC(t) = C(t) - C(t-), and © implies that 



dC(t) 



(7) 



0, for t^t k ; 
R k , for i = t k . 

The effectiveness of Di Paola and Falsone's formula expressed by © has 
been verified extensively [161 El HI 12 • 



3 Derivation of the alternative formula 



Define Y(X, t) as 



J! 



(8) 



i=i 

where , t) (j = 1,2, •••) are defined by ©. Using the notation Y 

defined in (IH1), we can write the Stratonovich SDE ([5]) as 



dZ(t) = f(Z(t),t) dt + Y(dC(t),t). 
Construct an ordinary differential equation(ODE) as 

g{Z(t-) + y{X),t) for 0<\<dC(t), 



0) 



d»(A) 



(10) 



dA 

y(o) = o. 

We claim that Y(dC(t),t) as defined in © is the solution of the above 
ODE ([TO]) at A = dC(t). The proof is given as follows. 

First, we check the existence and uniqueness of solution for (|10p . It 
is well known that if the function g(x, t) is Lipschitz continuous, then the 
existence and uniqueness of (I10|) can be guaranteed. Since © requires the 
function g(x,t) be analytic with respect to x [3H], which implies Lipschitz 
continuity, the existence and uniqueness of the solution follows. 

Take Taylor expansion of g(Z(t— )+y(A, t),t) at A = 0, then (Tl0j) becomes 

dy(X) 



dA 



00 Qk 
k=l 



A=0 



A^ 



(11) 
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where 

z(X) = Z(t-)+y(X). (1 
In the following, we show by induction that 

—g( z (X),t)=g^(z(X),t), for k>l. (] 
For k = 1, we obtain from (1120 and (jlOh that 

^(AM) = (^(A)) (^(A),t)) 

A^O (^ WA),<)) 

^(A),t))(^| X y^(A),t)) 

(2) (^(A),t)), (1 



= 5 

which means that (|13p is true for fe = 1. 
Suppose (jX3j) is true for k = n, i.e., 

SA™ 

then we get 

Qn+1 



g(z(X),t)=g^ +1 \z(X),t)), (] 



5 (z(A),t) = ^(J^(z(A,t)) 
= |s ( " +1) WA),t)) 



= ^(A),t))(^| X y<7 (n+1) WA),t))) 

= ^+ 2 )(z(A),t)), (] 

which means that (|13p is also true for A; = n + 1. Now the proof of (|13l) 
complete. 

Substituting (fT3l) into (fTT]> . and using the fact that 

^ +1 )(z(A),t))| A=0 = <?(" +1 )(Z(t-),t)), (] 



we get 



dy(A) 



oo 



X k-1 



dX 



(k-l)\ 



(18) 



k=l 



Integrating from to A at both sides of (|18|) . we get 



oo 



A* 
k\ 



(19) 



k=l 



Comparing ()19p with ([8]), we see the equivalence 



Y(dC(t),t)=y(dC(t)) 



(20) 



where the right-hand side y(dC(t)) is the solution of the initial value problem 
pop at A = dC(t). In other words, we obtain an equivalent and alternative 
expression for the Stratonovich SDE © or (|9|) 



where y(dC(t)) is the solution of the ODE (pH at A = dC(t). 

4 Some remark and examples 

As shown in section 3, when g(x,t) is analytic, Y(dC(t),t) defined in ([8]) is 
equivalent to y(dC(t)) obtained from the ODE (fTUj) . i.e., the two forms of 
SDEs, ([5]) and (f2T|) are equivalent. However, compared with ©, (f2~Tj) has 
significant advantages in the following aspects. 

First, for a given g(x,t), one immediate advantage of the alternative 
expression f)21 1) over ([5]) is that it is much easier to obtain analytical solution 
from the ODE (|10p than finding the sum of the infinite Taylor series in ([5]). 
Second, the proposed alternative expression (12ip is applicable under much 
more general condition than the original one (J5J) . For example, when g(x,t) 
is only continuous or piecewise continuous rather than analytic, (|2ip is still 
valid while ([5]) fails to exist. More importantly, the new expression (I21D 
could have significantly improved performance in numerical implementation 



dZ(t) = f(Z(t),t)dt + y(dC(t)) 



(21) 



6 



than the original expression ([5]). The form in (JSj) is actually the Taylor 
expansion of the solution of (|10|) . which requires computing the high-order 
derivatives of the function g(x, t) and possibly take a large number of terms 
to achieve high accuracy. However, since the alternative form f)21 1) is given 
in form of ODE, it can be solved by numerous algorithms such as Runge- 
Kutta or multistep methods, which could have much better convergence 
and stability properties than Taylor expansion methods. The readers are 
referred to any numerical analysis book such as [5] for more details about the 
Runge-Kutta, multistep, or Taylor expansion methods for solving ordinary 
differential equations numerically. 

We note that the alternative expression, as given in (|2ip . turns out to be 
a special case of the Marcus integrals, which have been defined for general 
semimartingale with jumps. For more details about the Marcus integral and 
its application, readers are referred to [TO l ITT | |9| IT]. 

Example 

In this example, f(x,t) and g(x,t)x are chosen such that the analytical 
solution of the SDE is available to compare with the numerical results. 
Take f(x,t) = x, and g(x,t) = x, the SDE ([5]) becomes 



where Y(dC(t)) is expressed, by Di Paola and Falsone's formula, as 



or by the alternative expression, as the solution of the following ODE at 



dZ(t) = Z(t)dt + Y{dC{t)) 



(22) 




(23) 



k=l 



X = dC(t), 



dy(A) 



Z(t-) + y{\) for 0<X<dC{t) 



(24) 



< 



2/(0) = 0. 



V 



It is easy to check that the theoretical solution of (|22p is e 



t+C{t) 
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In the simulation, the compound Poisson process C{t) is taken as 

N{t) 

C(t) = J2 R kU(t-t k ). (25) 

k=l 

wher N(t) is a Poisson process with intensity parameter as 10 and R k (k = 
1, ■ ■ ■ ,N(t)) are normally distributed random numbers independent of each 
other. 

When no jumps occur, Y(dC(t)) = and We solve (j22j) by second order 
Runge-Kutta method with time step as 0.01. We compute the jump of the 
solution Z(t) by Di Paola and Falsone's formula and the alternative formula, 
respectively. For the former, we truncate the infinite consequence by taking 
the first 6 terms, and for the latter, we solve ODE (I24p numerically by 
uisng second order Runge-Kuta method with step size equal to 0.1. Figure 

1 shows one random path of the Poisson compound process, and Figure 

2 compares the corresponding computational results with the theoretical 
one. It can be seen from Fig. 2 that when the jumps are computed by the 
alternative formula, the numerical result agrees very well with the theoretical 
one. However, when using Di Paola and Falsone's formula, the numerical 
result deviates from the theoretical one significantly, especially when the 
jump size is relatively large. We note that the accuracy for Di Paola and 
Falsone's formula could be improved by keeping more terms when truncating 
the infinite series. However, the number of terms needed to achieve the 
desired accuracy with Di Paola and Falsone's formula is not easy to predict. 
The alternative formula is preferred partly because its accuracy can be easily 
controlled by the step size and the order the Runge-Kutta or other methods. 
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Figure 1: One random path of the driving process 
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Figure 2: Comparison of the numerical results with the theoretical one 
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